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Abstract. In a recent article (Gances, Deleurence and Lewin, Commun. 
Math. Phys. 281 (2008), pp. 129-177), we have rigorously derived, by means 
of bulk limit arguments, a new variational model to describe the electronic 
ground state of insulating or semiconducting crystals in the presence of local 
defects. In this so-called reduced Hartree-Fock model, the ground state elec- 
tronic density matrix is decomposed as 7 = 'y'^^^ + Qv,sp, where -fp^j. is the 
ground state density matrix of the host crystal and Qiy,i!p the modification 
of the electronic density matrix generated by a modification u of the nuclear 
charge of the host crystal, the Fermi level sp being kept fixed. The purpose of 
the present article is twofold. First, we study more in details the mathemati- 
cal properties of the density matrix Qi/,ep (which is known to be a self-adjoint 
Hilbert-Schmidt operator on L^(R'^)). We show in particular that if J^^ u 0, 
Qv,ep is not trace-class. Moreover, the associated density of charge is not in 
I/^(R^) if the crystal exhibits anisotropic dielectric properties. These results 
are obtained by analyzing, for a small defect v, the linear and nonlinear terms 
of the resolvent expansion of Qi/,ep. Second, we show that, after an appropri- 
ate rescaling, the potential generated by the microscopic total charge (nuclear 
plus electronic contributions) of the crystal in the presence of the defect, con- 
verges to a homogenized electrostatic potential solution to a Poisson equation 
involving the macroscopic dielectric permittivity of the crystal. This provides 
an alternative (and rigorous) derivation of the Adler- Wiser formula. 



1. Introduction 

The electronic structure of crystals with local defects has been the topic of a 
huge number of articles and monographs in the Physics literature. On the other 
hand, the mathematical foundations of the corresponding models still are largely 
unexplored. 

In [3], we have introduced a variational framework allowing for a rigorous char- 
acterization of the electronic ground state of insulating or semi-conducting crystals 
with local defects, within the reduced Hartree-Fock setting. Recall that the re- 
duced Hartree-Fock (rHF) model - also called Hartree model - [53] is a nonlinear 
approximation of the A'^-body Schrodinger theory, where the state of the electrons 
is described by a density matrix 7, i.e. a self-adjoint operator acting on L'^{E?) (in 
order to simplify the notation, the spin variable will be omitted in the whole paper) 
such that < 7 < 1 and whose trace equals the total number of electrons in the sys- 
tem. The rHF model may be obtained from the usual Hartree-Fock model (TB] by 
neglecting the so-called exchange term. It may also be obtained from the extended 
Kohn-Sham model [6] by setting to zero the exchange-correlation functional. 

Our variational model is derived from the supercell approach commonly used 
in numerical simulations, by letting the size of the supercell go to infinity. It is 
found [5] that the density matrix of the perfect crystal converges in the limit to a 
periodic density matrix ^'^^^ describing the infinitely many electrons of the periodic 
crystal (Fermi sea). In the presence of a defect modelled by a nuclear density of 
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charge v, the density matrix of the electrons converges in the limit to a state 7 
which can be decomposed as [3] 

(1) 7 = 7por + Qi'.e-p 

where Qv,eF accounts for the modification of the electronic density matrix induced 
by a modification v of the nuclear charge of the crystal. The operator Qu.e^ depends 
on V as well as on the Fermi level ep, which controls the total charge of the defect. 

Loosely speaking, the operator Qu,ef describing the modification of the electronic 
density matrix should be small when u itself is small. But it is a priori not clear 
for which norm this really makes sense. The mathematical difficulties of such a 
model lay in the fact that the one-body density matrices ^^^^ and 7 are infinite 
rank operators (usually orthonormal projectors), representing Hartree-Fock states 
with infinitely many interacting electrons. 

Following previous results on a QED model [TUl HH H] describing relativistic 
electrons interacting with the self-consistent Dirac sea, we have characterized in [5] 
the solution Qv,eF © the minimizer of a certain energy functional E^^^-^ on 
a convex set K, that will be defined later. This procedure leads to the information 
that 

(2) e 62, Q-,- , Q++ e Si. 

In the whole paper we use the notation 

Q 7pcrQ7pcr Q ^ 7p(,j(3(l — 7pcr) 

Q^- (1 - 7p"or)Q7p°or (1 - 7p°or)Q(l " 7pV), 

and we denote by 61 and S2 respectively the spaces of trace-class and Hilbert- 
Schmidt operators on L^(R^). A definition of these spaces is recalled at the begin- 
ning of Section [2.31 for the reader's convenience. 

Property ^ implies that the operator Qu,ew is compact, but it does not mean a 
priori that it is trace-class. This mathematical difficulty complicates the definition 
of the density of charge. For Q G 61, the density of charge can be defined by 
Pq{x) = Q{x,x) where Q{x,x') is the integral kernel of Q; it satisfies Pq = 
Tr((5). Let us emphasize that these formulae only make sense when Q is trace- 
class. 

In [3] we have been able to define the density pq^ associated with Qv^e^ by a 
duality argument, but this only gave us the following information: 

(3) PQ^e e i (R ) and / / — dxdx < 00. 

"'^"^ Jr^ iR3 |a; — x I 

Also, following [TUl and using ^ one can define the electronic charge of the state 
counted relatively to the Fermi sea j^^^. via 

Tro(Q.,,,) := Tr(Q++ ) + Tr(Q--J. 

It can be shown that when is small enough (in an appropriate sense precised 
below) , 

Tro(Q.^ep) = 0, 

hence the Fermi sea stays overall neutral in the presence of a small defect. 
In [3], we left open two very natural questions: 

(1) is Qv,eF trace-class? 

(2) if not, is pg^ nevertheless an integrablc function? 

The purpose in the present article is twofold. First, we prove that Qv^sf '^^ ricvcv 
trace-class when Jj^a u 0, and that, in general, pq^ is not an integrable function 
(at least for anisotropic dielectric crystals). These unusual mathematical properties 
are in fact directly related to the dielectric properties of the host crystal. They show 
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in particular that the approach of ^ involving the complicated variational set /C 
cannot a priori be simplified by replacing /C with a simpler variational set (a subset 
of Si for instance). 

In a second part, we show that our variational model allows to recover the Adler- 
Wiser formula [TJ |^ for the electronic contribution to the macroscopic dielectric 
constant of the perfect crystal, by means of a homogenization argument. More 
precisely, we rescale a fixed density ly G L^(M'^) n _L^(R'^) as follows 

iy^{x) := rfv{r^x), 

meaning that we submit the Fermi sea to a modification of the external potential 
which is very spread out in space. We consider the (appropriately rescaled) total 
electrostatic potential 



(77 ^x) 



of the nonlinear system consisting of the density v^j and the self-consistent variation 
PQ^^ ep of the density of the Fermi sea. We prove that WJ converges weakly to Wi, , 
the unique solution in iS'(M^) of the elliptic equation 



(4) -div (emVI^^) = 47ri/ 



where Em is the so-called macroscopic dielectric permittivitij^, a 3 x 3 symmetric, 
coercive, matrix which only depends on the perfect crystal, and can be computed 
from the Bloch-Floquet decomposition of the mean-field Hamiltonian. As we will 
explain in details, the occurence of the dielectric permittivity Em, or more precisely 
the fact that in general Em 7^ 1, is indeed related to the properties that Qv,eF is not 
trace-class and pg^ is not in L^(]R'^). 

This article is organized as follows. In Section [21 we briefly present the reduced 
Hartree-Fock model for molecular systems with finite number of electrons, for per- 
fect crystals and for crystals with local defects. In Section [3] we study the linear 
response of the perfect crystal to a variation of the effective potential, the nonlin- 
ear response being the matter of Section 16.31 Note that the results contained in 
this section can be applied to the linear model (non-interaction electrons), to the 
reduced Hartree-Fock model, as well as to the Kohn-Sham LDA model. We then 
focus in Section [4] on the response of the reduced Hartree-Fock ground state of the 
crystal to a small modification of the external potential generated by a modification 
ly of the nuclear charge. We prove that for 1/ small enough and such that v 0, 
one has Tro(Q,y,£p) — while the Fourier transform pq^ ^^{k) of Pq^ does not 
converge to when k goes to zero, yielding Qu,ef ^ ®i- We also prove that if 
the host crystal exhibits anisotropic dielectric properties, pQ^ (fc) does not have 
a limit at fc = 0, which implies that pq^ ^ L^(R'^). Finally, it is shown in Sec- 
tion [5] that, after rescahng, the potential generated by the microscopic total charge 
(nuclear plus electronic contributions) of the crystal in the presence of the defect, 
converges to a homogenized electrostatic potential solution to the Poisson equation 
Q involving the macroscopic dielectric permittivity of the crystal. All the proofs 
are gathered in Section [6] 



2. The reduced Hartree-Fock model for molecules and crystals 

In this section, we briefiy recall the reduced Hartree-Fock model for finite sys- 
tems, perfect crystals and crystals with a localized defect. 



To be precise, it is only the electronic part of the macroscopic dielectric permittivity, as we 
do not take into account here the contribution originating from the relaxation of the nuclei of the 
lattice (the nuclei are fixed in our approach). 
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2.1. Finite system. Let us first consider a molecular system containing Af non- 
relativistic quantum electrons and a set of nuclei having a density of charge p""'^. If 
for instance the system contains AI nuclei of charges zi , • • • , € N \ {0} located 
at - ,Rk& R^, then 

K 

fc=i 

where mi, - • • ^uik are probability measures on R^. Point-like nuclei correspond 
to mfc — 6 (the Dirac measure) while smeared nuclei are modeled by smooth, 
nonnegative, radial, compactly supported functions ruk such that J^^ nik = 1- 

The electronic energy of the system of Af electrons in the reduced Hartree-Fock 
model reads [131 13j 

(5) fp„.c(7) = Tr (-^A7) - J^^ p, (p™^ * | • T^) + \d {p„ p,) . 

The above energy is written in atomic units, i.e. fi,= l,m=l,e = l and = 1 
where m is the mass of the electron, e the elementary charge, h the reduced Planck 
constant and Eq the dielectric permittivity of the vacuum. The first term in the 
right-hand side of ([5]) is the kinetic energy of the electrons and !?(•,•) is the classical 
Coulomb interaction, which reads for / and g in Ls(R'^) as 

(6) / MiM...,.4,/ 2«i*, 

Jr3 Jr3 \x - y\ iR3 \kY 

where / denotes the Fourier transform of /. Here and in the sequel, we use the 
normalization convention consisting in defining /(fc) as 

/(fc) = (2^)-i / /(:r)e-*-dx. 
JR? 

In this mean-field model, the state of the M electrons is described by the one-body 
density matrix 7, which is an element of the following class 

T'-^ = I7 e 5(i2(R3)) I < 7 < 1, Tr(7) = Tr (V^7%/^) < ooj, 

S[l?{^y) denoting the space of bounded self-adjoint operators on L^(R'^). Also 
we define Tr(— A7) := Tr(\/~A7-\/— A) which makes sense when 7 g . The set 
is the closed convex hull of the set of orthogonal projectors of rank M acting 
on L^(R'^) and having a finite kinetic energy. 

The function p~^ appearing in ([5]) is the density associated with the operator 7, 
defined by p-y{x) = ^{x, x) where 7(0;, y) is the kernel of the trace class operator 7. 
Notice that for all 7 € 7-"^, one has > and .^p^ e H^{M?), hence the last two 
terms of © are well-defined, since p^ e ^^(R^) n ^^{M?) C hiim?). 

It can be proved (see the appendix of [23j) that li Af < X^fcli (neutral or 
positively charged systems), the variational problem 

(7) inf{£p„uc(7), 7eP^} 

has a minimizer 7 and that all the minimizers share the same density p^. 

2.2. The perfect crystal. The above model describes a finite system of A/" elec- 
trons in the electrostatic field created by the density p""'^. Our goal is to describe 
an infinite crystalline material obtained in the bulk limit Af 00. In fact we shall 
consider two such systems. The first one is the periodic crystal obtained when, in 
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the bulk limit, the nuclear density approaches the periodic nuclear distribution of 
the perfect crystal: 

fQ\ „nuc nuc 

Ppgjr being a 7?.-periodic distribution, TZ denoting a periodic lattice of M.^ . The 
second system is the previous crystal in presence of a local defect: 

(9) r^^pii'i + u. 

The density matrix 7p(,j. of the perfect crystal obtained in the bulk limit ([5]) is 
unique It is the unique solution to the self-consistent equation 



(10) 



(11) 



/per -^( — ooieplv pcr7 



where V^or is a 7?.-periodic function satisfying 

-Al/per = 47r {pI^^ - pI"^) , with pI^^{x) = 7°„(a;, x), 

and where ep G R is the Fermi level. The potential Vpor is defined up to an additive 
constant; if V^cr is replaced with V^or + C, ep has to be replaced with ep + C, in such 
a way that ^^^j. remains unchanged. The function V^cr being in Lpgj.(R'^), it defines 
a A-bounded operator on i^(R^) with relative bound zero (see [TH Thm XIII.96]) 
and therefore H^^^ is self-adjoint on L?(R^) with domain H'^{E?). Besides, the 
spectrum of iJpgr is purely absolutely continuous, composed of bands as stated in 
PI Thm 1-2] and [11 Thm XIII. 100]. 

More precisely, denoting by TZ* the reciprocal lattice, by F the unit cell, and by 
F* the Brillouin zone, we have 

where for all q € F*, {£n,q)n>i is the non-decreasing sequence formed by the eigen- 
values (counted with their multiplicities) of the operator 

acting on 

Ll^,{T) := {u e iL(K') I u 7e-periodic} , 
endowed with the inner product 

{u,v)l2^^ ^ J uv. 

We denote by (un.?)n>i an orthonormal basis of Lp^^.{T) consisting of associated 
eigenfunctions. The spectral decomposition of {Hp^^.)q thus reads 

oo 

(12) (i7°^,)q = ^e„,5|u„,g)(u„,,|. 

Recall that according to the Bloch-Floquet theory [TF, any function / G L^(R^) 
can be decomposed as 

.fix)^£f,{x)e^''-^dq, 
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where f^, is a notation for |r*| ^ Jp, and where the functions fq are defined by 

(13) = + R)e~''-^''^"^ = ^ E + ^y"'^- 

Ren ' ' Kew 

For almost aU q e M^ fq £ ^^^.(r). Besides, fq+K{x) = fq{x)e-'^-'^ for aU K eU* 
and almost all g £ M'^. Lastly, 

ll/lli2(R3) = £^ ll/«lli^„(r) dq. 

If the crystal possesses N electrons per unit cell, the Fermi level ep is chosen to 
ensure the correct charge per unit cell: 

(14) N = J2\{qer* \sn^q<ep}\. 

n>l 

In the rest of the paper we will assume that the system is an insulator (or a semi- 
conductor) in the sense that the N*"^ band is strictly below the {N + 1)*^' band: 

:= maxe^v.g < min e^v+i.^ := ^n+i- 

In this case, one can choose for ep any number in the range (E^,E^_,_j^ ). For 
simplicity we will take in the following 

= 2 

and denote by 

9 = ^N+i - > 

the band gap. 

2.3. The perturbed crystal. Before turning to the model for the crystal with a 
defect which was introduced in [5], let us recall that a bounded linear operator Q 
on L^(K^) is said to be trace-class [HI [l^ if J2i {fi^ VQ*Qfi) 1^2 < 00 for some 
orthonormal basis {ipi) of L^(R'^). Then Tt{Q) = iVijQ'Pi) is well-defined 
and does not depend on the chosen basis. If Q is not trace-class, it may happen 
that the series {^i,Qfi) 1^2 converges for one specific basis but not for another 
one. This will be the case for our operators Qu,eF- 

A compact operator Q = J2i'^i\'Pi){fi\ ^ iS(L^(M'^)) is trace-class when its 
eigenvalues are summable, |Ai| < 00. Then the density 

+00 

Pq{x) = Q{x,x) = YX^\ipt{x)\'^ 

i=l 

is a function of L^(R'^) and 

-I- 00 

PQ- 



+00 „ 
Tr(Q) = ^A,=== / 



On the other hand, a Hilbert-Schmidt operator Q is by definition such that Q*Q is 
trace-class. 

We now describe the results of ^Sj dealing with the perturbed crystal. We have 
proved in [5] by means of bulk limit arguments that the ground state density matrix 
of the crystal with nuclear charge density Pp"Jr + ly reads 

7 = 7por + Q^'.EF 
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where Qu,eF obtained by minimizing the foUowing energy functional 

(15) £;,,,,(Q)=Tr(|i/0„-eF|(Q++-0— ))- / Pq(^^ * I • P') + ^^(pq, Pq) 
on the convex set 

(16) /C = {OeQ| -7°er<Q<l-7pV} 
where 

(17) Q = {Qe62|0* = Q, O^^eSi, Q++e©i, 

|V|Q e 62, |V|Q— |V| e 61, |V|Q++|V| e 61}. 

Recall that 61 and 62 respectively denote the spaces of trace-class and Hilbert- 
Schmidt operators on L^(M'^) and that 

Q 7pcrQ7pcr Q .— 7p(,^(5(l — 7pcr) 

Q + ~ := (1 - 7°er)Q7°er (1 " 7p°or)Q(l ^ 7p^)- 

It is proved in [3j that although a generic operator Q e Q is not trace-class, it 
can be associated a generalized trace Tro(Q) = Tr((3++) + Tr(Q^^) and a density 
PQ G L'^{M.^) n C where the so-called Coulomb space C is defined as 

C:=|/e5'(M3)|i?(/,/)<oo where i^(/, /) 47r ^^^^ rffcj . 

Endowed with its natural inner product 

C is a Hilbert space. Its dual space is 

C := {V e L^{R^) I WV e (i2(M3))3} ^ 
endowed with the inner product 

{Vi,V2)c' ■■= ^ [ VVi-VV2^^f \kfvAk)V2{k)dk. 

Note that if Q e /C n ©1, then of course Tro(Q) = Tr(Q), pq{-) = Q{-, •) e L^{M.^) 
and Tr(Q) = pq. 

The energy functional i?i/.ep is well-defined on JC for all v such that (i^* | • |^^) € 
L'^{S?) +C'. The first term of -Ejy^ep makes sense as it holds 

Ci(l-A) < |ijO,,-£F| <Ci(l-A) 

for some constants < ci < Ci < oo (see [SJ Lemma 1]). The last two terms of 
are also well defined since pg € L^(IR3) n C for all Q e /C. 
The following result is a straightforward extension of Theorem 2 in 3J , allowing 
in particular to account for point-like nuclar charges: if i/ is a Dirac mass, • S 

L2(M3) 

Theorem 1 (Existence of a minimizer for perturbed crystal). Let v such that 
{v-k\- \~^) £ L^(IR3) -\- C' . Then, the minimization problem 

(18) inf{£;,,,p(Q), Qg/C} 

has d TniTiiuiizGT Qi/^s^, and all the uiiTiiinizeTs of 111 8\) shuve the sauie density p^y ^p. 
In addition, Qu.e^ '^^ solution to the self- consistent equation 

(19) Q.,,p = l(-oo,ep) {Hi,, + (p.,ep -^)*\- - l(-oo,ep] (i^pcr) + ^ 

where S is a finite-rank self-adjoint operator on ^^(M^) such that < ^ < 1 md 
Ran{5) C Ker {H^,, - ef) ■ 
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Remark 1. Our notation Qu,ep does not mean that minimizers o/iJ^ ^p are nec- 
essarily uniquely defined (although the minimizing density p^^^F itself unique). 
However, as we will see below in Lemma\^ when v <^ \ in an appropriate sense, 
one has 5 — Q hence Qi,^ew indeed unique. 

In this approach, the electronic charge of the defect is controUed by the Fermi 
level ep, and not via a direct constraint on Tro(Q) (see [3] for results in the latter 
case). When 

the minimizer Qu^s^ uniquely defined. It is both a Hilbert-Schmidt operator 
(Tr((5^_jp)^ < oo) and the difference of two orthogonal projectors (since (5 = 
then). In this case, Tro(Q,y,£p) is always an integer, as proved in [TUl Lemma 2]. 
The integer Tro(Qi',eF) can be interpreted as the electronic charge of the state 
7 = 7per + Qv,eF (mcasurcd with respect to the Fermi sea 7por)- We will see later 
in Lemma [5] that Tro((3i/,ep) = whenever v is small enough. 

Note that the fact that Qi/,ep is both a Hilbert-Schmidt operator and the differ- 
ence of two orthogonal projectors automatically implies that the generalized trace 
of Qv,ei? is well-defined since 

with Qjjp > and Q^^p < 0. Let us remark incidently that the condition 
Tr((5^gp) < c» is required in the Shale-Stinespring Theorem which guarantees 
the equivalence of the Fock space representations i4j defined by ^^^^ and 7 = 
7por + Qy,eF respectively. 

One of the purposes of this article is to study in more details the operator Q^.e^ 
and the function p^^e^ ■ 

3. Linear response to an effective potential 

In this section, we study the linear response of the electronic ground state of a 
crystal to a small effective potential V G L'^{M.^) +C'. This means more precisely 
that we expand the formula 

Qv = l(-oo,eF] (^pcr + V) - 1(-oo,£f] (^per) i 

in powers of V (for V small enough) and state some important properties of the 
first order term. The higher order terms will be studied with more details later in 
Lemma m Obviously, the first order term will play a decisive role in the study of 
the properties of nonlinear minimizers. 

As mentioned in the introduction, the results of this section can be used not 
only for the reduced Hartree-Fock model considered in the paper, but also for the 
linear model and for the Kohn-Sham LDA framework. In the reduced Hartree- 
Fock model, the effective potential is y = {pi^.e-p ~ ^) * I ■ 1^^- In the linear model, 
the interaction between electrons is neglected and V coincides with the external 
potential: V — Vcxt = —v-k \ ■ |^^. In the Kohn-Sham LDA model, 

V = (pj;.?p^ --)*!• r + -L°^P°e. + pJ;.?.^) - v^FHpU 

where v^c is the LDA exchange-correlation potential and p]f^^ the variation of the 
electronic density induced by the external potential Vcxt = —iy-k\ ■ \~^, see [4]. 

Expanding (formally) Qy in powers of V and using the resolvent formula leads 
to considering the following operator 

(20) Qi,v = ^ £ - ^Por)"' V{z- i7p",,)"' dz, 
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where is a smooth curve in the complex plane enclosing the whole spectrum of 
iJpg^ below ep, crossing the real line at ep and at some c < inf (T(_ffp(,j.). In order 
to relate our work to the Physics literature, we start by defining the independent 
particle polarizability operator xo- 

Proposition 1 (Independent particle polarizability). If V E L^(M'^) + C , the 
operator Qi,v defined above in (j20p is in Q and Ttq(Qi,v) — 0- If V E i^(R^), 
Qi,v G Si and Ti{Qiy) = 0. 

The independent particle polarizability operator xq defined by 



XoV := PQ, ^ 



is a continuous linear application from L^(S.^) to L^{M.^) and from L'^{1SP) + C to 
L2(K3) nc. 

The proof of Proposition [1] is provided below in Section 

For the cases we have to deal with, we can consider that the effective potential 
V is the Colomb potential generated by a charge distribution p: 

V ~ p-k\ ■ for some p. 

We will have p ~ —v for a linear model (non-interacting electrons) and p = p^.ep —v 
for the nonlinear reduced Hartree-Fock model. Following usual Physics notation, 
we denote by the Coulomb operator: 

■vc{p) = I ■ r\ 

which defines an isometry from C onto C. If p £ v-'^{I?i^-^)+C') = v^^ {I? {;^)) ^- 
C, then we have Vc{p) € L'^(M.^)+C', hence Qi.v,{p) e Q and Pq,^^^^^-^ e L'^{R^)r\C. 
We now define the linear response operator 



(21) C{p) := -pQ,, 



and we concentrate on the study of the operator C. As £ = — Xo^c it follows from 
Proposition[l]that £ maps C into C ni^(M'^). The reason why we have put a minus 
sign is very simple: in the rHF nonlinear case, we will have 

where r2 contains the higher order terms, and and we will rewrite the above equality 
as 

(22) {1 + C){iy~ p,,,^)^u-h. 

This motivates the following result. 

Proposition 2 (Self-adjointness of the operator C). C defines a bounded nonneg- 
ative self-adjoint operator on C. Hence 1 + C, considered as an operator on C , is 
invertible and bicontinuous from C to C. 

The latter property will be used in Section [H 

We have considered the linear response for all reasonable V^s (or p's). We now 
assume that V ^ p-k\-\^^ with a density p £ L^{M.^) and we derive some additional 
properties of C{p). Note that as L'^{R^) C v-'^{L'^(M.^)) + C, we have v^p) G 
L^(M'^) +C'. The following statement is central in the mathematical analysis of the 
dielectric response of crystals. 

Proposition 3 (Properties of C{p) when p e L^). Let p e _L^(M'^). Then, C{p) £ 
L^(IR'') nC, C{p) is continuous on M.^ \ TZ* , and for all a £ S"^ (the unit sphere of 



(23) hm /:(p)(77a) = (a^ia)p(O) 
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where L eM?^^ is the non-negative symmetric matrix defined by 

2 

a N +00 

(24) VfceM^, k^Lk^^S2 Y: 



■ ^9, 



where the £n,q 's and the Un,q 's are the eigenvalues and the eigenvectors arising in 
the spectral decomposition of {H^^^)q. 
Additionally, 

(25) Lo = iTr(L) > 0. 

Proposition [3] shows that C{p) is not in general a function of L^(]R'^) even when 
p G L^(KP), as when L ^ Lq (i.e. when L is not proportional to the identity matrix), 

C{p) is not continuous at zero (note that L ~ Lq characterizes isotropic dielectric 
materials). However the following holds: for any radial function ^ G C^(R^) such 
that 0<^<1, ^=lon 5(0, 1) and ^ = on \ b{0, 2), we have 

(26) lim [ C{p){x)£_{R-^x)dx = Lo [ p. 

4. Application to the reduced Hartree-Fock model for perturbed 

crystals 

Let us now come back to the reduced Hartree-Fock framework and the decay 
properties of minimizers. Our main result is the following 

Theorem 2 (Properties of the nonlinear rHF ground state for perturbed crystals) . 
Let V G L-^(R'^) n L^(R'^) he such that J^^ v ^ and | • I^^Hl^+c small 
enough. Then the operator Q^^e^ satisfies Tro((32y,ep) — hut it is not trace-class. 
If additionally the map L : ^ defined in Proposition [21 is not constant, then 
Pi,,eF >^o^ in L^{M!^). 

The proof of Theorem [2] is a simple consequence of our results on the operator 
C stated in the last section, and of the continuity properties of higher order terms 
for an density p. A detailed proof is provided in Section [^751 

As previously mentioned, the situation L = Lq characterizes isotropic dielectric 
materials; it occurs in particular when 7?. is a cubic lattice and Pp^Jr has the sym- 
metry of the cube. For anisotropic dielectric materials, L is not proportional to the 
identity matrix, and consequently p,y,£p ^ L"'^(M'^). 

Formula (pS)) for Lq is well-known in the Physics literature [TJ [5S]. However 
to our knowledge it was never mentioned that the fact that Lq > is linked to 
the odd mathematical property that the operator Qi,^ew trace-class when 

/r3 ^ ^- The interpretation is the following: if a defect with nuclear charge v 
is inserted in the crystal, the Fermi sea reacts to the modification of the external 
potential. Although it stays formally neutral (Tro((5iy,ep ) = 0) when v is small, 
the modification p^.ep of the electronic density generated by 1/ is not an integrable 
function such that J^^ Piy,eF — 0, as soon as J^^ v ^ 0. 

For isotropic dielectric materials, L = Lq, and we conjecture that the density 
Pi/^ep is in L^(R'^). In this case, one can define the total charge of the defect 
(including the self-consistent polarization of the Fermi sea) as /jja {v — p^.ew ) ■ For 
v small enough, the Fermi sea formally stays neutral (Tro((5iy,ep) — 0), but it 
nevertheless screens partially the charge defect in such a way that the total observed 
charge gets multiplied by a factor (1 -t- Lo)~^ < 1: 



1 + Lo' 
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This is very much similar to what takes place in the mean- field approximation of no- 
photon QED [121 IS]- In the latter setting, the Dirac sea screens any external charge, 
leading to charge renormalization. Contrarily to the Fermi sea of periodic crystals, 
the QED free vacuum is not only isotropic (the corresponding L is proportional 
to the identity matrix) but also homogeneous (the corresponding operator C has a 
simple expression in the Fourier representation) , and the mathematical analysis can 
be pushed further: Gravejat, Lewin and Sere indeed proved in ^ that the observed 
electronic density in QED (the corresponding py.ej,) actually belongs to i-^(R'^). 
Extending these results to the case of isotropic crystals seems to be a challenging 
task. 

When L is not proportional to the identity matrix (anisotropic dielectric crys- 
tals), it is not possible to define the observed charge of the defect as the integral 
of v Pv.sy since Pu,eY not an integrable function. Understanding the regularity 
properties of the Fourier transform of p^^e^ is then a very interesting problem. In 
the next section, we consider a certain limit related to homogenization in which 
only the first order term plays a role and for which the limit can be analyzed in 
details. 

5. Macroscopic dielectric permittivity 
In this section, we focus on the electrostatic potential 

(27) V = {y-p,^,,)^\-r^ 

generated by the total charge of the defect and we study it in a certain limit. 
We note that the self-consistent equation can be rewritten as 

(28) u-p,^,^ = {l + C)-^v-{l + C)-^f2- 

Therefore for the nonlinear rHF model, the linear response at the level of the density 
is given by the operator (1 + C)^^ . We recall from Proposition [5] that £ > on C 
and that (1 -|- £)~^ is a bounded operator from C to C. 

In Physics, one is often interested in the dielectric permittivity which is the 
inverse of the linear response at the level of the electrostatic potential, i.e. 

e-^ := v,[l + C)-\-^ . 

Note that can be recast into 

(29) V ^e-\,{v)-v,{l+Ly^h- 

A simple calculation gives (we recall that xo is the polarizability defined in Propo- 
sition (TJ which is such that C ~ — Xoi'c) 

Therefore one gets 

(30) e-^ = l + v,{l + C)-\o- 
We also have 

(31) e^v,{l + C)v-^ 
which yields to the usual formula 

(32) £ = 1 - VcXo- 

The basic mathematical properties of the dielectric operator e are stated in the 
following 
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Proposition 4 (Dielectric operator). The dielectric operator e = 1 — VcXo is an 
invertible bounded self-adjoint operator on C , with inverse e"^ = 1 + Uc(l + /^)^^Xo- 
The hermitian dielectric operator e — Vc ^ evi is an invertible bounded self- 
adjoint operator on L^(R'^). 

The proof of Proposition is a simple consequence of the properties of xo and C, 
as explained in Section [6.91 

Even when v e i^(R'^), applying the operator (1 + C)^^ creates some discon- 
tinuities in the Fourier domain for the corresponding first order term (1 + £)~^i/ 
in Equation (pS)) . If we knew that the higher order term r2 was better behaved, 
it would be possible to deduce the exact regularity of p^.e^- We will now con- 
sider a certain limit of (j28p by means of a homogenization argument, for which the 
higher order term vanishes. This will give an illustration of the expected proper- 
ties of the density in Fourier space at the origin. For this purpose, we fix some 
ly e i^(R^) n L^(R^) and introduce for all 77 > the rescaled density 

v,j{x) := rfv{T]x). 

We then denote by VJ' the total potential generated by i^^, i.e. 

(33) K'':=K-P.,„ep)*|-r\ 
and define the rescaled potential 

(34) W:f{x):^r^-'v:j{r'x)- 

Note that the scaling parameters have been chosen in such a way that in the absence 
of dielectric response (i.e. for e^^ = 1, 7~2 = 0), one has WJ] — Vc{i^) — v*] ■\'^^ for 
all r)>0. 

Theorem 3 (Macroscopic Dielectric Permittivity). There exists a 3 x 3 symmetric 
matrix em ^ 1 such that for all v £ L^(R'^) n L^(]R''), the rescaled potential 
defined by ^34^ converges to W^, weakly in C when 77 goes to zero, where Wv is the 
unique solution in C' to the elliptic equation 

-div (emVWi.) = 47r7^. 

The matrix Em is proportional to the identity matrix if the host crystal has the 
symmetry of the cube. 

From a physical viewpoint, the matrix is the electronic contribution to the 
macroscopic dielectric tensor of the host crystal. Note the other contribution, orig- 
inating from the displacements of the nuclei [18j . is not taken into account in our 
study. 

The matrix £m can be computed from the Bloch-Floquet decomposition of i?per 

as follows. The operator = Vc^^^s~^vl^^ commuting with the translations of 
the lattice, i.e. with for all i? G 7?., it can be represented by the Bloch matrices 

ii^K^K' i<l)]K,K'eTZ' )qer' ■ 

V/eL2(R3)^ f^f^q + K)= e-'^,{q)f{q + K') 

K'eK' 

for almost all q E T* and K E TZ* . We will show later in Lemma [5] that eK.K'iv^) 
has a limit when 77 goes to 0+ for all fixed a E S'^. Indeed one has 

lim £0 oirjcr) = 1 + La 

7;-»0+ 

where L is the 3x3 non-negative symmetric matrix defined in . When K, K' 7^ 
0, £k,K'{'>1^) has a limit at 77 = 0, which is independent of a and which we simply 
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denote as Sk,k'{0)- When K = Q but K' ^ 0, the limit is a linear function of a: 
for all K' en*\ {0}, 

lim eo,K'{w) = Pk' ■ cr, 

for some Pk' G C^. The electronic contribution to the macroscopic dielectric per- 
mittivity is the 3x3 symmetric tensor defined as f5] 

IfcP 

(35) VfceM^ k^eyik^ lim , , , , . 

By the Schur complement formula, one has 
1 



£o,K{'rik)[C{rik) ]k,k'£k' ,o{vk) 



where C{rjk)~^ is the inverse of the matrix C{r]k) = [sKK'{vk)]K.K'eTZ'\{o}- This 
leads to 

lim ,^ = \k\^ + k^Lk- V (/3K-fc)[C(0)-i]^f,K'(/3i^'-fc) 

where C(0)^^ is the inverse of the matrix C(0) = [sKK'{0)]K.K'en'\{o}^ hence to 



(36) 



K,K'en'\{o} 



As already noticed in [5], it holds 

1 < £m < 1 + i. 

Formula (|36p has been used in numerical simulations for estimating the macro- 
scopic dielectric permittivity of real insulators and semiconductors Uni HH [71 [5] ■ 
Direct methods for evaluating £m, bypassing the inversion of the matrix C(0), have 
also been proposed [5D1 [TS] . 

6. Proofs 

In this last section, we gather the proofs of all the results of this paper. 
6.1. Preliminaries. Let us first recall some useful results established in [3]. 

Lemma 1 (Some technical estimates). Let A be a compact subset ofC \ a{Hp^^). 

(1) The operator B{z) :— {z — Hp^^.)^^(l — A) and its inverse are bounded 
uniformly on A. 

(2) The operators |V| x \z~Hp^j.\^2 and iVKz — H^^^.)^^ are bounded uniformly 
on A. 

(3) There exists two positive constants < ci < Ci < oo such that 

(37) ci(l~A)<|i/°„-£F|<Ci(l-A). 

(4) IfVe L2(M3) +c', hpa,.,V] e ©2 and there exists a constant C G R+ 
independent of V such that 

(38) l|[7pcr,^]l|e. <q|^||L^(K3)+C'. 

Besides, if V e L'^(R^) for some 1 < g < oo and ifWe L'P{R^) for some 
I < p < oo, then 

\Her,V]\\^^<C\\VVh.iM^y 

We denote as usual by &p the space of all operators A such that Tr(|A|P) < oo, 
endowed with the norm ||A||ep := Tr(|74|'')p. 
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Proof. For (1), (2), (3) and the first assertion of (4), see the proofs of [3, Lemma 1 
and O Lemma 3]. The last estimate is obtained like in [31 p. 148] by writing 




It then suffices to use the Kato-Seiler-Simon inequality (see [5T] and [211 Thm 4.1]) 

(39) Vp > 2, ||/HV).9(a;)||e, < (2^)"^||5||lp(e3)||/|Up(k3) 

and the fact that B(z) is uniformly bounded on □ 

6.2. Proof of Theorem [H Let f be such that V = {ly * \ ■ \-^) e L'^{m.^) + C . 
As C^iJS?) is dense in L^(M'^) and is included in C, V can be decomposed for all 

> Q &sV ^V2,^ + with V2,r, e L2(R3), V^^ G C and 11^2,^^^ < r/. Denoting 
by v'^ = -(47r)"^AV"j^, we obtain v'^ e C and 

VQeQ, -/ pQV>-rj\\pQ\\L2-D{:.'^,pQ). 

By [21 Prop. 1], we know that there exists a constant C £ M+ such that 

vgeQ, \\pq\\l^-<C\\Q\\q. 
Besides, for all Q e /C, < Q++ - Q^^ with g++ > and Q-- < 0. Hence, 

VQ e /C, ||pq|U2 < C + C'Tr ((1 - A)(Q++ - Q— )) . 
Using (|37p and choosing 77 > such that 277C" < ci leads to 

(40) VQe/C, i?,,,,(Q)>|-Tr((l-A)(0++-Q— ))-C'ry-il?(;.;,i/;). 

The above inequality provides the bounds on the minimization sequences of (jl8p 
which allow one to complete the proof of Theorem [1] by transposing the arguments 
used in the proof of [3, Theorem 2] . □ 

6.3. Expanding Qy- In this section, we explain in details how to expand 

Qv l(-oo,eF] {Hpci- + V) - l(-oo,eF] (^pcr) ) 

and give the properties of each term in the expansion. The multiplicative operator 
associated with some V € L^(R'^)+C' is a compact perturbation of H^^j., so that the 
operator H^^^. + V is self-adjoint on i^(M''). When V is small, it is then possible to 
expand Qv in a perturbative series, using the resolvent formula. For this purpose, 
we consider a smooth curve in the complex plane enclosing the whole spectrum 
of iJper below ep, crossing the real line at ep and at some c < inf (T{Hp^^). We 
furthermore assume that 

d{a{H°J,A) = ^ where A = {z € C | ^) < |} , 

d denoting the Euclidian distance in the complex plane and g the band gap (see 
Fig.lH). 

The following result will be useful to expand Qv- 



Note there is a sign misprint in the corresponding formula at the top of p. 148 in 
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A 



Figure 1 . Graphical representation of a contour C C enclosing 
a{Hp^^.) n (—00, ep] and of the compact set A. 



Lemma 2. There exists a > such that if V E L^(M^) + C is such that 

||V^|U2+c' < a, 

then 

(41) a + y) n A - 0. 

Moreover ep i a + V), Qv e Q and Tto{Qv) = 0. 

Besides, there exists an orthonormal basis {ip^)i>i of the occupied space = 
i?a77.(7pj,j.) and an orthonormal basis [tf'l)i>i of the occupied space TL+ ~ Ker{'-^^i^^) 
such that in the orthonormal basis {{^^), {^t)) of L^(R^) — Tl-(B'H+, 

diag{-ai, -02, ■ ■■) 11 diag{bi,b2, ■ ■ ■] 



(42) Qv 
with 



diag{bi,b2,- ■ ■) 11 diag{ai,a2,- ■ ■) 



+ CXD 

E 

=1 



< fli < 1, \^ a,; < 00, bi — y/ ai{l — Oi). 



The meaning of ([1T|) is the following. As mentioned above, any F e L^( 
C defines a compact perturbation of H^^j.; hence the essential spectrum of the 
Hamiltonian remains unchanged: 

(Tess (-ffpcr + V) = (Toss (-ffpcr) • 

This in particular means that only eigenvalues of finite multiplicity may appear in 
the gap (S^, E^_|_j^), and they can only accumulate at or For V small 

enough in +C', these eigenvalues will be localized at the edges of the gap, i.e. in 
a vicinity of and It can be seen that the charge TioiQtv) jumps as t is 

increased when an eigenvalue crosses the curve "^lo and that it is a constant integer 
when this does not happen. By continuity, we deduce that Tto{Qv) = 0: for V 
small enough, no electron-hole pair is created from the Fermi sea. 

The representation (|42p of Qv was proved in [llj and it can be interpreted in 

terms of Bogoliubov states. Each 2x2 submatrix ( ] can be seen as a 

\ bi ai J 

mrtea/ electron- hole pair. A reaZ electron- hole pair would be observed for Oi — 1. It 

1 X 

is easy to see that the eigenvalues of Qy (including multiphcities) are {—a^ , )i>i. 
Thus a necessary and sufficient condition for Qy being trace class reads J^'^i < 00. 
We now provide the 

Proof of Lemma\^ Let C G M+ be such that for aU 2; e A, ||i?(z)|| < C (see the 
first statement of Lemma [T]) . For all z e A, 

z - -V={z- H%,){1 - B{z){l - A)-iy). 
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As L'^{R^)+C' C L'^{R^) + L^{R^), it follows from the Kato-Seiler-Simon inequality 
([55)) that there exists a constant C such that 

Vy e l2(m3)+c', ||(1-A)-i^|| <C'\\V\\l2+c'- 



As V defines a compact perturbation 



ii\\vh2+c' < (cc')-\a(i/pV. + y)nA 

of -ffpor: it follows from a standard continuity arguments that for ||V^||l2+c' < 
(CC")~^, the set a{Hp^^. + V)r] (—00, ep] lays inside the contour yielding 



Thus, 



1 
2i7r 



pcry 



Besides, still in the case when || V^Hl^+c < {CC 



H 



per 



so that 



< 



m c^c'\\v\\l2 



2n l^CC'WVh^+c' 

We now set a = ((1 + {27r)-^ C\'^\)CC')-\ For all V G L^{M.^) + C such that 
||F||z,2+c' < a, it holds <t{H°^^ + n A = and ||(3v|| < 1. We conclude using 
[TUl Lemma 2] and [H Theorem 5]. □ 



We are now going to expand Qv using the resolvent formula. We already know 
that 



7d 



1 



1 



(-00, ep] 



2i7r jtjf 

It now follows from the proof of Lemma [2] that 



dz. 



1 



1 



for all e £2 



C such that ||T^||l2^c' < Qf; yielding for such V^'s 



1 



2nr 7^ 



dz. 



One important result of the present section is the following 

C such that II 1^11 L2+C' < o;- 



Lemma 3 (Resolvent expansion). Let V G ^^(M'' 
T/ien, for allK gN\{0}, 



(43) 
where 



(44) 
(45) 



Qy — Qi, 



1 



Q 



K,V 



!K+1,V 



V{Z- ifp"er) 



fe 



dz 



Qky = —(piz^H^-V) '\v{z-H^^^ 



dz 



For all k > 1, the operator Qk.v in Q and Tro(Qi;.v') = 0. For all k > 1, the 
operator Qky is in Q and Tro(Qfe,y) — 0. For all fc > 6, the operators Qk,v (^i^d 
Qk,v o.f^ trace-class and Tr((3fc,y) = 0. 
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Note that by linearity, the operators Qk,v are well-defined for all V € L'^{M!^)+C' , 
and not only for small Vs. It can in fact be shown using the same arguments as 
in the proof of Lemma [3] that for all fc > 1, 

{Vu--- ,Vk)^^(f {z- (z - H^,,)-' ■■■Vk{z- H^,,)-' dz 

J ^ 

is a continuous fc-linear application from (L^(M'^) to Q. 

Let us now detail the 

Proof of Lemma It follows from the proof of Lemma [5] that each term of the 
expansion makes sense in the space of bounded operators on L^(R'^). We now 
have to prove that Qk,v and Qk,v are in Q and that their generalized trace is equal 
to zero. We start by noticing that Qv is indeed a minimizer for the functional 

E{Q) Tr(|i/0,, - £f|(Q++ - Q—) - [ Vpg 

on /C. Theorem [T] with the nonlinear term erased then implies that Qv G Q- 

Let us consider Qiy. Decomposing V as V = V2 + V with V2 € L^(R^) and 
V e C C i'5(R^), and using the Kato-Seiler-Simon inequality ([Ml) and the first 
assertion of Lemma [1] Qi.Va ^ ®2 and Qiy S Gq. Hence, Qiy is well-defined 
in 6q. a straightforward application of the residuum formula then shows that 

HIT Jc^ 

we can make use of Lemma [T] to conclude that Qiyr E &2- Obviously, the same 
holds true for Q^y,, so that Qiy, and henceforth Qiy, are in 62- As |V|(z — 
Hpcr)~^ is a bounded operator, uniformly in 2: € ^, it is easy to check that | VlQi^y^ 
and |V|Qi,v' both are Hilbert-Schmidt operators. Finally, |V|Qi,y G 62 and 
therefore Qiy € Q. As Qiy = Qiy = 0, we obviously get Tro((3i,y) = 0. 

Let us now consider Qky for fc > 2. The potential V being in L^(M^) + L^(M^), 
we have 

• Q2y e 63 and |V|Q2,y e 63; 

• Qky € 62 and |V|(5fe,y € S2 for all fc > 3. 

As usual [iniE]) the next step consists in introducing 7°^,^ + (7por)^ = 1 in (^1]) in 
places where (-ffpor ^ •^)^^ appears, and in expanding everything. We will use the 
notation 

and similar definitions for all the other terms. A simple application of the residuum 
formula tells us that Q^y'^ = Q^y^ — 0. Therefore, Q^y = Q2v^ and 1/ = 
Q2v^- Now we remark that the terms Q^y~ and Q2v^ involve two terms of 
the form Tper^'^CTpor)^ = bpcn ^Klpcr)''" (o^ adjoint) in their formula. Using 
Lemma HJ we obtain that Q2V, Qtv^ \'^\Q2v\^\ ^^"^ I^IQ^i/l^l are trace-class 
operators. Likewise, Q^y, y, iVlQ^TylV] and |V|(5^y|V| are trace-class oper- 
ators. Lastly, Q^y = Q^y^ + Q^y^, both operators of the right-hand side involv- 
ing one term of the form 7per^(7per)^ = [ipei-^VKlper)^ ■ Consequently Q^^ and 
jVlQj^y are Hilbert-Schmidt. Repeating the same argument for Q^y, we obtain 
that Q2y and |V|(52,y are Hilbert-Schmidt. Therefore, all the operators Qky are 
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in Q. As Qv also is in Q, Qk,v G Q for all k > 3. It also follows from the first 
assertion of Lemnia[T]and the Kato-Seiler-Simon inequality that Qk,v is trace-class 
for A; > 6. 

Using ([55]) and the Kato-Seiler-Simon inequality ([55]) . we then easily obtain that 
for all fc > 1, there exists a constant Ck S K+ such that 

yV e L\R^)+C', WQkMlQ < Ck\\V\\l.+c' 

and that for all fc > 2, there exists a constant Ck G K+ such that 

yV e L'{R^)+C' s.t. \\V\\l2+c' < «, IIQfe^ylls < CkWVWl^+c- 

Let V e L'^(R^)+C' such that ||V^||l2+c' < a. For aU t e [0, 1], UtT^IU^+c' < " and 

Qtv = Qi.ty + • • • + Qkav + QK,tv = iQi,v + • • • + t^QK,tv + QK+i,ty- 
As we know that Tro(Qty) — 0, we obtain that for all t E [0, 1], 

= <Tro(Qi^y) + • • • + t''TTo{QK,tv) + T^roiQx+uv) 

with |Tro(Q/^+i,ty)| <j\QK+iMQ < Cii+it'-'+^VW^+^c,. Hence, TvoiQk^v) - 
for aU fc > 1 and Ti-Q{Qky) = for aU fc > 2. □ 

6.4. Proof of Proposition [TJ Let V G L^{W^) + C. We already know from 
Lemmainithat Qi.y G Q and that Tro(Qi,y) = 0. Decomposing V as V = V2 + V 
with V2 G L-^(M'^) and V £ C\ and proceeding as in Section [Ql we obtain 

\\Ql.V.\\Q<Cm\L-^, 

\\Qi,v'\\Q<C\\V'\\a. 

We infer that xo is a continuous linear application from from L^(R'^)+C' to L^(M'^)n 
C. 

Let us now examine the case when V S L^(R'^). Using again the Kato-Seiler- 
Simon inequality and the first assertion of Lemma[Tl we obtain G 61 and 

||Qi,y||ei < C\\V\\li. 

Consequently, xo defines a continuous linear application from L^(M'^) to i^(R'^). 
As from the residuum formula, Qiy = Qiv = 0, we get Ti{Qiy) =0. □ 

6.5. Proof of Proposition m As p e ^^(R^) C v^'^{L'^{R^)) +C, we have S 
L^(R'^) n C. The operator C can be explicitely calculated in Bloch transform. We 
start from the Bloch-Floquet decomposition of -ffp^r- fo'" / G iJ^(R'^), 

{H"^,j){x) ^ £aKM)e"''' dq 

where (see Eq. 

+00 

(-^pcr)<3 ~ y ] £n,q\Un,q) {Un,q\- 
n=l 

Note that by time-reversal symmetry, 

'^n,—q — '^n^qj ^71, ~q — ^n^q- 

Denoting by U = Vcip), we can write the Bloch matrix of the operator Qiy as: 

(46) [Q^,v]qq, ^^(f{z- {H^,,),rWq^q,{z - (ff^ej,')"' dz . 
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Inserting the spectral decomposition of H^^^ in we obtain 



per 

N +00 



(47) [Ql,v]qg' ^ E ( "I {Un,q,Vq-q'Un',g') ^2 (r) l^n.g) | 

„=1„'=JV+1 

^n' ^q ^n.q' 

Remark 2. /n the following we will write series of the form (j47p and we will 
invert sums and integrals without giving any justification. To see that such a series 
is absolutely convergent, one can use the fact that there exists a and b in Rj^ such 
that for all n> I and q GT* , 

£n,q > an^^'^ — b. 

This bound is easily obtained by comparison with the eigenvalues of the periodic 
Laplacian. It follows that there exists C G M+ such that 

C 

\{Un,q,Vq^q'Un' ,q') j^2^^^y,^ \ < -Jyy^ 

for alll <n< N, all n' > N + 1 and all q, q' e T* . 

If an operator A d Q has a Bloch matrix Aqq> , then we have 

(48) {pA)q{x) ^ £ Aq>^q,^q{x,x)dq'. 

This formula is obtained by writing, for any real- valued function / S L^(]R'^), 
{iPA)qJq)L2^^^j^^dq = j pA{x)f{x)dx 



Tr{Af) = j^TvLi^^(r){{Af)q,,q,)dq' 
dq j- dq" Tri2^^(r) (A,/, 



dxj- dq' j- dq Aq,+q^q>{x,x)fq{x). 



We deduce that 



(49) [C{pMx) = f dq'Y, 

( ^^n.q' ^ {,P ~^ \ ' I ^q^^i-' ^q' — q) 1^2 fp-\^n,g' ('^)'^n',ij'— q(*^) 

\£n' ,q' — q £n,q' P^'^ 

~t~ {^n'.qS (P ^ I ' I ^q'^n^q' — q^ 1^2 (Y)^n' .q' {x)Un^q' —q(^X^ | . 

^n' ,q' £n,q' — q ^'^^ J 

The next step consists in decomposing the operator L, as the sum of a singular part 
and a regular part, corresponding respectively to the low and high Fourier modes 
of the Coulombic interaction kernel | • |^^. More precisely, we choose some smooth 
function ^ which equals 1 in a small neighborhood i?(0, S) of and outside the 
ball B(0, 2(5), with (5 > such that 5(0, 25) C V* . Then we define 

where J-~^ is the inverse Fourier tranform. Similarly we define 
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and note that Vr E L^{M.^). This being said, we have by the Kato-Seiler-Simon 
inequaUty that 

hence Cr{p) G L^(R'^) and /^a Cr{p) = by Proposition [1] Consequently, Cr{p) € 

C°(R3) and CA?)(0) = 0. 

Let us now deal with the singular part of C{p)- Using the definition of the 
Bloch-Floquet transform, we obtain that 

and that for almost all g G F*, 



{T-\^p))^{x) = 
This implies that for almost all q G F*, 



(27r)? 



(50) 



{p*Vs)q{.x) = 47r 



(27r)t £,{q)p{q) 



Therefore we get for almost all g G F*, 
(51) 
where 



C,[p)^{x)^{2^)i^p{<l) 



N +00 



(52) 



I I "'r* n=ln'=N+l 



'TL'.q'—q ^n,q' 



('^n,g' : ^n' .q' —q} 1^2 ^p^'^n,g'(*^)'^n',ij'— q (-^) 



,q ^n.q' ~q 

It follows that almost everywhere in F 



^ {'^n' m' : '^7l,q' —q} T 2 (y\^n'q'(,'^^'^n,q'—q{-^^ ]• 

£n',q' ~ £n,q'~q P"^ ' J 



Cs{p){q) = (2^)-5 j Cs{p)q{x)dx 



where 
(53) 



B{q) 



N +OC 
n=l n'=N+l 



dq' 



{'^n,q' T "^n' ,q' ~q) j-p-) 



We now remark that the above formula may be written 



Stt 



Biq) = ^—Mq)TTL2^^ 



'€ 



dz -f dq 



(7° ) h° )^ 



q'^q 



r* ^ (-^pciOg' ^ (-^pcr)ij'-g 



We recall [17] that g 1— > (7por)^ is a smooth periodic function and that (7per)^ is 
for all q a rank-A^ orthogonal projector. It is then easy to deduce that q ^ B{q) 
is a continuous periodic function on R"^. Consequently, Cs{p) and therefore C{p) 
are continuous on F* \ {0}. Using similar arguments, one obtains that C{p) is 
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continuous on R^\TZ*. In order to study the limit of Cs{p){q) when q goes to zero, 
we use the relation 



(54) 



V|2 



{Un,q',Un'q'-q)Ll^^(J') 

to rewrite jCs{p){w) for S and r] > small enough as 
where 



N +00 



n=l n'=N+l 



{a ■ VUn,q', Un' ,q' -no) Ll^ 



(r) 



dq' 

r* {^n' ,q' —r}cr ^n,q'^{^n' ,q' —rja ^n,q' "I" VQ^ ' ^ 2 ) 



(recall that ^ = 1 in the vicinity of 0). Again the above formula may be rewritten 



as 



1^ I n=l 



n,q' 



which shows that when 77 goes to zero, L^{(j) converges to a^La while p{r]a) con- 
verges to /o(0). 

We now turn to the proof that Lq > 0. We note first that 



^o = 3pE E 

I I n=l n'=N+l-"' 



(Vu„,5',«„',g')i2 (D 



i^^n' ,q' ^n^q'^^ 



hence Lq = would imply {'Vun,q,Un',q) L^^^ir) = for all q G F*, all n < N and 
all n' > + 1. Hence for i — 1,2,3, dx^ would stabilize the space Xg spanned by 
{ui^q, UN,q) for any q. Next we differentiate the eigenvalue equation for Un,q and 
get 



A 

-— ~iq-V + — + Vper - en,q 



From this we deduce that (Sxj^er) would also stabilize Xq. This means that we 
would have 



( '^\,qWl 



\ UN,q{x) 



UN,q{x) 



for some matrix Aq depending only on q. As u\fi{x) > for all x (it is the first 
eigenfunction of a Schrodingcr operator), wc deduce that {dxiVpci){x) would be for 
any a; G F an eigenvalue of Aq. By continuity and periodicity we infer that V^er 
would be constant. This is in contradiction with the assumption that the host 
crystal is an insulator or a semiconductor. □ 
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6.6. Proof of Proposition [2l The proof of Proposition [T] shows that £ defines a 
bounded linear operator on C. Besides, for all pi and p2 in C, 

mpi),P2) - / [/:(pi)](x)K(P2)](x)dx 

1 



= (p Tr {{z - ffO„)-i«e(pi)(z - i?p",,)-ii;c(p2)) dz 

= mp2),pi) = {piXiP2))- 

Therefore, C is self-adjoint on C. Lastly, denoting hy V ^ ^cip), we have for all 

pec, 

(£(p),p)c= / C{p)ix)V{x)dx^ j- I C{p),{x)Vq{x)dxdq 




-n.q 

2 



> 0. 



We conclude that £ is a bounded positive self-adjoint operator on C. Consequently, 
1 + £ : C ^ C is invertible. □ 

6.7. Expanding the density to higher orders. In the previous sections, we 
have studied the first order density —£{p). For the proof of our Theorem [2] on the 
reduced-Hartree-Fock model, we need to consider the higher order terms. Each of 
the operators Qk,v (for fc > 1) and Qk.v (for k > 2) defined in Lemma [5] being in 
Q, the expansion (j43p can be rewritten in terms of the associated densities, yielding 
the following equation in L^(]R^) n C: 

(55) PQv = PQi.v H ^ PQk.v + PQk+i,v ■ 

We introduce the quadratic operator r2 defined by 

which is continuous from v'^^{L'^{E?)) + C to i^(M^) n C, and the nonlinear map 
fa from v-^{Bo,) to L'^{W^) n C {B^ denoting the baU of L'^{M.^)+C' of radius a), 
defined by 

We obtain 

(56) r{p) := pq„^,^, = -C{p) + r^ip) + f3(p). 

The next lemma is concerned with the second and third order terms of the expansion 
([561). We will assume that p G ^^(R^) n L^{9?). Then p G Li{W^) C C, so 
that Vc{p) G C and |l'yc(p)|lc' < C'IIpIIlirl^, where C is a universal constant. In 
particular, there exists a constant 7 > 0, such that 

WpWl^hl^ < 7 \W{p)\\l^+c' < \\vc{p)\\c' < a. 

Lemma 4 (Nonlinear terms in the expansion). Let p G L^{R^) n L^(]R'^). Then 

(1) Q2{p) £ 6p for all p > 1 and the Fourier transform r2{p) of r2{p) is 
continuous on and vanishes at k = ; 

(2) If in addition, \\p\\L^nL^ < 7; then Qsip) G 61, r3{p) G L^(M'^) and 

Tr(Q3(p))- / h{p)^0- 
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Proof of Lemma^ As p E (R^) O L'^ (M.^) , wc deduce from Young inequality that 
y = I . |-i is in LP{R^) for 3 < p < c» and that is in (^'(K^))^ for aU 
I < g < 6. Therefore, 

for all I < g < 6, by Lemma [T] Arguing as in the proof of Lemma [31 we obtain 
that Q2,v & &p for all p > 1. 

We now concentrate on the regularity of r2 (p) — pq^ v at the origin. Arguing 
like in the proof of Proposition [31 we only have to study the density associated with 
the operator 

R2:^^([){z- H^,,r' ip^v,) (z - i/p",,)-i {p*v,){z- H'pJ-'dz. 

We will for simplicity only treat the term i?^^"*", the other ones being similar. 
Following the proof of Proposition [31 we obtain 

N 



(57) p^(g) = (2^)-i^ i '^^'i 

n=lm>N + lm'>N+l''^' 



>N + 1 

{'^n,q' 7 (P * ^s^q' —r'^m.r} {^m.r i (P * ^s^r — q'+q^m'.q' — q} {'^m' ,q' —q : ^n,g' ) 
{,^71, q' ^m,r)(£"n,g' ^m' .q' —q') 

Changing r ^ r — q' and using as before (|5D)) . we see that for \q\ small enough, 
(58) p^+{q)=AV2iY. E J2 i dq'-f drC(-r)C(r + q)x 

^ p{ r)p(r4-^) {Uji^q' 7 '^rn^r^q' ) {'^m,r-\-q' 5 ^m' .q' ~q) {'^m' ,q' —q •> '^n.q'} 
|r|2|r + q\^{en,q' - £m,r-g')(^n.9' " £m',q'-q) 

Next, using (|5l)) . we obtain 

n=l m>Ar+l „i'>JV+l ■'i * •'i* III ^1 

|p(-r)| |(Vlt„,q', 



|£m,r+g' - Eri.g' - r • (r + 2g')/2| 

|p(r + g)| VM„,g')| 



{^7n.r—q' ^n,q'^{^m' ,q' —q ^n,q'^\^n,q' ^m',q'—q ~^ Q ' {,Q '^Q )/'^l 

Note that choosing the support of ^ small enough we have 

\e„7,r+q' - e„,g' - r • (r + 2g')/2| > c(m2/3 + 1) 
uniformly for r, (7' G F* and n — I..N . Similarly, taking q small enough we get 

\en,q' - em',q'-q + q ' [q - 2q')/2\ > c((m')'/^ + 1). 
Using these estimates and the fact that p E we deduce that 

and the result follows. 

To establish that Qsy is trace-class, and therefore that r3(p) is integrable, we 
write 

Qsy = Q3,v + Q4,v 

and proceed as above to prove that both operators in the right hand side are trace- 
class. As TioiQs.v) — 0, we readily conclude that Ti'iQs.v) — /rs ^^(p) =0. □ 
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6.8. Proof of Theorem [2j We now have all the material for proving Theoreni[51 

The first step is to confirm that if the external potential v-k\ ■ is small, so is 
the effective potential [v — Pi^^sf) * I • I hence the results of Section [6731 can be 
applied. 

Lemma 5. There exists (3 > such that if 

ik*i-r'iu=+c' <A 

then 

where a is the constant used in the formulation of Lemma [H Consequently, the 
solution to U8\) is unique and satisfies Tro(Q,/,ep) = and 



per 

1 



dz. 



(60) = 7r-(P (2-i?p^-(p.,ep-i^)*|-r') -(^-i?p°er) 

Proof of Lemma[5i Let < S < 1 and ly such that ||!^* | • < S. This 

implies that 

with V2 e L^iR^), v' e C, \\V2\\l^ < S and \y\\c < S. We then deduce from (jiOl) 
that 

= £;.,ep(0) > i?.,e.(g.,ep) > yTr ((1 + |Vn(Q++ - Q-- )) - C"5- ^. 
It follows that there exists a constant c G independent of S and such that 

Tr((l + |V|2)(Q++ -Q--J) <c^. 
Using again the inequalities < — Q^^, Q^^ > 0, Q < 0, we obtain 

IIQ.,spIIq < 2c(5i 

Therefore, there exists a constant c' independent of 5 such that for all ly such that 
lk*|-|-^||L^+c'<'5, 

We obtain the desired result by choosing (3 — min 

The proof of Theorem [5] is a simple consequence of the results of Sections 16.31 
and 16. 71 We assume that | • I^^Hl^+c < /? in such a way that Lemma [5] can be 
applied. This gives us that Tro((5i/,£F) — and that — pi.^^^) *\ ■ |"^||l2+c' < 
Hence we can use the expansion of Lemma [5J 

If Qv,EF were trace-class, then we would have p„,j;p G L^(M'^) and 

Pv,eviO) = / Pi'^ep = Tr((3^.,ep) = 0. 



On the other hand, we would obtain from the expansion (|56p 

(61) Pi^.EF = -^{Pu..eF -v)+ r2{pu.,ej, - v) + r^iPu.ej, " 

By Proposition [3] and since we have assumed Pu.ev ^ i^(K^), we know that 
lim (p^{ria) + Cipy^e-e - i^)iw)] = -{cr"^ La)v{0) 

for all a £ S^. On the other hand we have by Lemma |4] that the Fourier transform 
of the second and third order terms r2(p,y,£p — i^) and r3(py_£p — v) vanish at the 
origin. It would then follow that {a"^ La)t'{0) — for all a G 5^, which obviously 
contradicts (|25p . Therefore, Q^^ep is not trace-class. 
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Let US know assume that Pu.e^ ^ -/^^(M'^). The same arguments lead to 

p,^ep(0) = -(<T^i<T)(p,^,,(0)-i7(0)) 

for aU (T G 5^. This is only possible if L = Lq. □ 

6.9. Proof of Proposition |4l Lete^l — UcXo- It follows from Proposition!!] that 
e is a bounded self-adjoint operator on C . Besides, using the fact that C :— — Xo^c 
we easily see that 

= 1 - v,xa + v,{l + Cy'xo + Veil + Cy'Cxo = 1- 

Likewise, e [1 + ■i;c(l + ^)'^Xo] = 1- Hence, e'^ = [l + Wc(l + Cy^Xo]^^ ■ 

Lastly, vi is an invertible bounded linear operator from L^(M'') onto C' . Besides, 
for ah / and g in L^(R^), 

{vc f^vig'j^^ = {f,g)L'^- 

As e is an invertible bounded self-adjoint operator on C , e is an invertible bounded 
self-adjoint operator on iy^(M^): 

{sL9)l^ = (ylef^vlg^^^ ^ (vlf.vieg^^^ = (/,?3>l2. 
The proof is complete. □ 

6.10. Proof of Theorem [3j Let v € L'^{^^) n L'^{^^). Introducing the dilation 
operator {Unf){x) — f{rix), we can write = rj^UnV. The operator [/^ is an 
isometry of L^(M'^) and satisfies J7^(p(/c) = ri~ ^ ip{r]^^ k) . It follows that 

I? (fc) = - 

yielding 

\WJc = V^\W\\c- 

Hence, for 77 > small enough, Iji^,, ★ | • |^^||l2^c/ < j3. Arguing like in the proof of 
Lemma [S] we obtain 

Therefore, 

We therefore may use the self-consistent equation ([25]) and get 

(62) Z^,, - Pix^.ep = (1 + Cy^l',^ - (1 + Cy'^hiVrj - Pi.„,ep), 

where r2 (p) : = r2 (p) -f- f 3 (p) . The bounds of the proof of Proposition [3] and the fact 
that (1 -|- is a bounded operator on C imply 

(63) \\{l + C)-y2{l^^-Pu,.e.)\\c<C\\v^-Pu,„eAl<C'lJ. 

For convenience, we are going to study equation ([5^ in L?{R^). We therefore 
introduce 

_i 1 
fr, ■■= Vc ^ WH and g := vi v. 

We note that is indeed bounded in L^(M'^) by the choice of the scaling in ■ 

Making use of the relation rjVcUrjVc ^ = Ujj, we can rewrite (l62|) as 

(64) /, = u; u^g - r^- u;vl{^ + cr^hi^n - P-,,e,), 
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1 _i 

where we recall that e^^ — Vc (1 + C)~^Vc ^ is a bounded self-adjoint operator on 
I/^(M'^). Our bound ([M]) on the nonlinear term shows that 

||/,-C/;e-i{/,.g|U2(R3)<Cr;i 

Hence the theorem will be proved if we show that U* e^^ C/,,g converges weakly in 
L^(R^) to the correct hmit as 77 ^ 0+. 

This will follow from the following two important lemmas. 

Lemma 6 (The macroscopic dielectric permittivity). We denote by iq : L'^^^{Y') 
Lp(,^(r) the Bloch transform of the operator i. We also denote by eo |r|^2 the 
(normalized) constant function and by Pq the orthogonal projection on {eo}^. The 
following hold: 

(1) The maps q ^ £q and q 1— > e"^ are continuous on T* \ {0} and uniformly 
bounded with respect to q. 

(2) For all a G 5^, £40-60 converges strongly in Lp^j.{T) to 6cr(a;), where for all 
k gM.^, the periodic function bk{x) is defined by 

(65) bk = {\k\^ + k^Lk)eo 

2^V^^i f ,,^f hlX' 



r(fc.V)t 



and where Gq is the operator defined on Lp^^{T) as 

^ ^ \/47r f]^ e^^^ ^ f e~^^^ 

^0/= J2 -tkt^i^ ""^^"^^ Ik= f{x)—^dx, 

Keiz'xm 1^1 1^1= |r|2 

and which satisfies PoG§ = G§Pq. 
(3) The family of operators PoiqPo seen as bounded self-adjoint operators acting 
on Poipgr(r) is continuous with respect to q and one has 

(Po£g-Po)|p^i2^^(r) ^ C 

strongly as q ^ 0, where C > 1 is the bounded operator on Po-^pcrC-T) 
defined by 

1 r ^ / (7" )^ 1 \ 

(66) Cf^f + 2Gif rfg^E LrfO . Un.q'G§f]u^ 

Jr* „=i \ ^^pcr7<?' / 

for all f G P^Ll^,{T). 



(4) One has for all a G S 
(67) lim (eo, £,7^60) 



1 



Using (j67p . we may now define the macroscopic dielectric permittivity as follows: 



k^EMk := + k^ Lk - (Po^fc, C'^Pobk). 



As Po^A: is linear in k, it follows that Em is a constant 3x3 symmetric matrix. 
Theorem [3] readily follows from 

Lemma 7 (Limit of the linear term). Let g be a fixed function in L'^{R^). Then 
U*e^^Urig weakly converges in L^(]R'^) as rj ^ 0^ to the function whose Fourier 
transform is given by 

m 



l + ^-(Po6^,C-iPofo^ 
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Assuming Lemma [51 we first write the 

Proof of Lemma^ As e^^ is a bounded operator on it suffices to show 

that 



lim IU:e-'U,s,s') = / it 

for two functions g,g' S L^(M'^) such that both g and g' have a compact support 
(say in a ball of radius R) . As the Fourier transforms of C/^ 5 and of Ujjg' have their 
support in the ball of radius Rrj, for rj small enough such that B(0,Rri) C F*, we 
have by the definition of the Bloch-Floquet transform 



(£-1(7^3, = / {{e-')keo,eo)U,,gik)Ur,g'{k)dk 



{{£ ^),fkeo,eo)g{k)g'{k) dk. 
The result then follows from (|67p and the dominated convergence Theorem. □ 

It now remains to write the 
Proof of Lemma O Using (|49p and time-reversal symmetry, we deduce that for all 

„ N +00 

(68) {ej){x)^f + 2f dq'J2 E 

(un',g'+q,Un,q'{Vc)q f) {Vc)q [u^Un' ,q' +q) (x), 

en'.,q'+q- en,q'\ ' ^l^A^) V ^ .yT^J/JV^ 

1 

where {vc)q is the convolution operator by the corresponding Bloch component, 
which just consists in multiplying the K^^ Fourier coefficient of a function by 
(47r) 2 |iir + The above formula can be rewritten as 



r ^ I 

J = f + 2{v,)lj-^^dq'Y,i-^ 



\ IpcrJq'+q 



(69) iqf = / + 2(^;,)| f dq' ^ { _l / 

\-"-pcr)q'+q ^ri,q' J 

We note that for any g G F* \ {0}, {vc)q is a bounded (indeed compact) operator on 

Lpg^(F) and that g 1— > {vc)q continuous from F* \{0} to C{L'^^^{T)). The continuity 
of Eq when q stays away from is therefore easy to verify. 
For / e Poi^er(r), we have 



where Gq is the operator on Lp^^{T) which multiplies the K*^ Fourier coefficient 
of a function / by y/inlq + K\~^ except the coefficient corresponding to X = 

which is replaced by zero. Note that Gq Gg in norm. Formula ((69)) then shows 
that C{q) :— Po£g-Fb|poL2 (r) is bounded and converges as g ^ to the operator 
C defined on Popper (r) by Obviously, C > 1 on P^L^ ^.(F). 
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Next, we have for q small enough 
{iq - l)eo 



^r- ^i„,r]^_^i \q\[£n'.q'+q- en,q') 



TjT^O + f dq }^ -TTT- ; {Un,q'Un',q'+q) 



where B{q) was defined before in ((53|) . Now we use (j54p and get 

TV +00 



{_^7i' ,q' -\-r]a ^n.q' ){^n' ,q' -{-rja ^7i.q' ' ^ 2 ) 

Hence as 77 ^ 0"^, er^creo converges strongly in L'^^j.(r) to b^- 

The last step is to use the Schur complement formula which tells us that 

(e-^eo eo) = - 

°' ° (£geo,eo) - (Fo£9eo,C((7)-iFoe<;eo) ' 

The above convergence properties yield 

lim, (So^eo, eo) = — — 577 ,p , „ , \ 

as was claimed. □ 
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